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The apparent lack of suitable astrophysical sources for cosmic rays with E > 10 19 - 7 eV (UHECRs) 
is the "GZK Paradox". We argue that whatever mechanism produces them must also account for 
events down to ~ 10 18 ' 7 eV, including their isotropy and spectral smoothness. This rules out galactic 
sources and distributed sources such as topological defects; Gamma Ray Bursts (GRBs) are unlikely. 
We are lead to identify the powerful radio galaxy Cen A, at 3.4 Mpc, as the probable source of most 
UHECRs observed at Earth today, and to estimate the extragalactic magnetic field to be ~ 0.3/iG. 
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Above an energy of about 10 197 eV, cosmic ray pro- 
tons suffer inelastic collisions due to photopion produc- 
tion from the cosmic microwave background W. Less 
than 20% of protons survive with an energy above 3 x 
10 20 (1 x 10 20 ) eV for a distance of 18 (60) Mpc §; 
ultra-high energy (UHE) nuclei and photons lose energy 
even more readily. Yet more than 20 cosmic rays (CRs) 
have been observed with nominal energies at or above 
10 20 ± 30% eV [| with the record Fly's Eye event hav- 
ing 3.2 x 10 20 eV [Q. There are no apparent suitable 
astrophysical sources within the GZK distance to which 
these events point nor is there any significant break 
in the spectrum in the GZK region, as would be expected 
if higher energy UHECRs are attenuated. 

The CR spectrum can be described by a series of power 
laws. At the "knee", E ~ 10 15 5 eV, the spectral index 
steepens from -2.7 to -3.0; at ~ 10 177 eV it steepens 
further to -3.3. At E ~ 10 185 eV the spectrum flattens to 
an index of -2.7 || and is consistent, within the statistical 
uncertainty of the data (which is large above 10 20 eV), 
with a simple extrapolation at that slope to the highest 
energies, possibly with a hint of a slight accumulation 
around 10 195 eV. For a review and references, see ||. 

The simplest interpretation of this data is that above 
E ~ 10 18 ' 5 eV a new population emerges which domi- 
nates the more steeply falling galactic population, and 
this new population has an approximately E~ 2 ' 7 spec- 
trum up to the highest observed energies. This inter- 
pretation is supported by AGASA's analysis of arrival 
directions ||: around 10 18 eV the angular distribution 
correlates with the galactic center and is consistent with 
a galactic origin, while at higher energy this anisotropy 
disappears. The smoothness of the observed spectrum 
above 10 18 5 eV suggests that a single mechanism is re- 
sponsible for these events. Otherwise, there is an ap- 
parently miraculous matching of spectra from different 
mechanisms, such that the total spectrum is smooth. 

We adopt here the following interpretation of the data: 
(i) A single population of events dominates the CR spec- 
trum from about 10 18 5 eV to the highest observed en- 
ergies. (It will be useful below to distinguish, within 
this population, between "low energy" , 10 18 7 — 10 19 5 
eV, and "high energy", E > 10 20 eV, UHECRs, avoiding 
transition regions.) This single mechanism proposition 
is extremely powerful. It means that the high statis- 
tics "low energy" UHECR data, which have previously 
been generally ignored, can be used to draw rather strong 
conclusions, (ii) The UHECR population has a smooth 
spectrum without a substantial discontinuity at 10 1 
eV. (iii) The UHECRs are nucleons, although alterna- 
tive models are mentioned and constrained as well. 

Above about 10 19 7 eV, protons display a rapidly 
falling attenuation length, plateauing at order 10 Mpc 
for energies above « 10 20 5 eV (see e.g. 0). Thus inde- 
pendently of the nature of the sources of UHECRs and 
of magnetic fields which may deflect or confine them, 
"low energy" UHECR protons with pathlengths up to a 
few Gpc contribute to the flux at Earth, whereas "high 



energy" UHECR protons only reach Earth if their path- 
length is of order few 10's Mpc or less. We denote the 
ratio of the average pathlengths for low and high energy 
UHECRs the accumulation factor, f acc . For uniformly 
distributed sources active over cosmological times, f acc ~ 
few Gpc/ few 10's Mpc w 100. If the UHECR sources 
produce a smooth spectrum, and if Earth is located in 
a "typical" environment, the observed spectrum should 
have an offset in normalization between low and high en- 
ergy UHECRs approximately equal to the accumulation 
factor, f acc ~ 100. 

There seem to be only three ways to account for 
face ~ 1 in the observed spectrum. 

(i) The contribution from sources within the GZK dis- 
tance (« 10 — 20 Mpc) and corresponding time (30 — 60 
Myr) is comparable to the contribution from all other 
sources active since redshift of order 1/2 added together. 
Galactic sources satisfy this condition trivially, however 
we will see that isotropy is difficult to reconcile with a 
galactic explanation. For extragalactic sources, this con- 
dition implies that our local sources are significantly more 
concentrated or more powerful than average. 

(ii) The spectrum at the source could have the inverse 
structure to that of the attenuation length. Some models 
appear to give a satisfactory accounting of the spectrum 
by this means, but only as a result of requiring a fit over a 
too-limited range of energy. For instance models in which 
the spectrum is governed by hadronization of quarks, 
e.g., the Z-burst, superheavy-relics, and topological de- 
fect models, have hard intrinsic spectra ~ E^ 1 . Thus 
over the limited energy range ~ iq 19 - 5-20 e y, where the 
attenuation length is rapidly dropping, the attenuated 
spectrum can be fit to the observed UHECR spectrum. 
However the spectrum of such models has the entirely 
wrong shape below the GZK energy, and at 10 18 7 eV, 
where the flux is well measured, the model flux is more 
than an order of magnitude too small. 

(iii) GZK energy degradation can be circumvented by in- 
voking new physics, but models involving GZK-evading 
messengers are constrained by this discussion: the nu- 
cleon spectrum at the source must be below that of the 
messenger particle by a factor w 1/ face, or else the accu- 
mulation of those protons generates an unobserved offset. 

A popular possibility considered recently is a com- 
bination of i) and ii), taking advantage of the overden- 
sity of matter in the local supercluster and hardening 
the spectrum [||. However the matter overdensity on 
the GZK scale falls far short of the f acc level required 
to remove the offset, and therefore these efforts fail to 
account for the spectrum below ~ 10 19 4 when the fit at 
higher energies is acceptable; Fig. 7 of |f| shows that the 
prediction is more than a factor of 2 too low at 10 19 eV. 

In ref. we showed that extragalactic magnetic fields 
can be of order a few tenths /zG rather than the much 
smaller values previously typically assumed in discus- 
sions of UHECR propatation, and we argued that the 
approximate isotropy of arrival directions could be due 
to even the highest energy CRs propagating diffusively. 
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Now we explore whether this is consistent with the three 
conditions outlined above: i) correct directional proper- 
ties, i.e., sufficient isotropy ii) a spectrum which is ap- 
proximately a smooth power law from ~ 10 18 7 — 10 20 5 
eV, and iii) flux of "low energy" UHECRs from the local 
source should be comparable to that of all other sources 
in the Universe. 

The nearby radio galaxy Centaurus A (NGC5128) 
proves to be an excellent candidate source. We show 
that the three conditions exclude virtually all other al- 
ternatives, leading us to conclude that Cen A is probably 
the source of most UHE cosmic rays observed on Earth 
today. We now summarize some results on diffusion in 
disordered magnetic fields and propagation in diffusive 
media which will be needed below. 

Magnetic fields are trapped in ionized matter whose 
turbulent flow leads to extensive restructuring of mag- 
netic fields over large and small distance scales. It is 
plausible to assume a Kolmogorov form for the mag- 
netic field power spectrum as a function of wavenumber, 
B(k) w Bk' 11 ^ 0@ in this regime; B denotes the 
magnitude of the field at scale A, the maximum scale of 
coherent correlations, generally expected to be ~ 0.1 — 10 
Mpc. The Larmor radius of a proton in a magnetic field 
B orthogonal to its motion is Rl = 110 Eiq/B^q, kpc. 

Diffusion has two qualitatively distinct regimes, de- 
pending on whether particles are trapped inside mag- 
netic subdomains or not, causing different functional de- 
pendence on energy of the diffusion coefficient D(E). 
Kolmogorov diffusion applies when <C A, i.e., E < 
I0 21 B^qXm P cgV ', so the power spectrum of the magnetic 
field irregularities is important. Using the Kolmogorov 
spectrum, jilt] obtained 
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As the particle energy increases, there is a transition to 
Bohm diffusion when Rl f» A and the diffusion coefficient 
is of order the Larmor radius times velocity. In order 
that either Bohm or Kolmogorov diffusion apply even for 
the Fly's Eye particle, we require B m gAm P c ^ 0.3. Nu- 
merical simulations [[L2| show that for Bohm diffusion, 
D w 3 R L c, i.e., D{E) « 0.1 E 20 /B fJ , G Mpc 2 /Myr. Since 
the prefactor in the diffusion coefficient is approximately 
the same for either Bohm or Kolmogorov diffusion at 
the highest energy of concern here, E20 ~ 3, we sim- 
plify the discussion below to only Kolmogorov diffusion; 
given the low statistics above 10 20 eV, this will be ad- 
equate. After infering the pathlength of propagation, 
we will check that it is long compared to the mean-free- 
path « D/c, verifying a posteriori the reasonableness 
of the diffusive approximation even for the highest en- 
ergy CRs. For the low energy UHECRs, D(10 18 7 eV) w 
1/8 Mpc 2 /Myr (A 2 lpc 0.3/B AlG ) 1 /3. 

When no CRs are produced by an isotropic source and 
propagate diffusively without energy loss, in an environ- 
ment with a spatially uniform diffusion coefficient D, 



their number density as a function of distance r from 
the source, and time t since emission, is 



n(r,t) 



"0 



in D i) 3 / 2 



-r 2 /[4D(E,B)t] 
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At a distance R from the source, the number density 
of CR's of a given energy increases rapidly with time, 
reaches a maximum we will call the diffusion front at 
t = R 2 /(6D), and then drops slowly with time ~ t~ 3 / 2 . 
Since D is an increasing function of energy, the diffusion 
front arrives earlier for higher energies. In order to avoid 
unobserved structure in the energy spectrum, we infer 
that the diffusion front of UHECR's from the dominant 
source must have already reached Earth, for all energies 

10 18.7-20.5 eV _ 

The number of particles with velocity c hitting a unit 
area in a unit time in a uniform gas of density n(r, i) is 
n(r,t)c. Due to the gradient in the number density with 
radial distance from the source, the downward flux at 
Earth per steradian as a function of angle 9 to the source 
is 



f(9,r,t) = (1 + a cos 9) 



n{r, f)c 



4tt 



where 



a 



D(E,B) 
f c dr 



2tc 
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If the source is collimated into back-to-back pencil 
beams, eqn (|2|) is approximately replaced by the super- 
position of two sources each with half the strength, sep- 
arated by a distance 2d w 2D/c along the direction of 
coUimation and with t — > t — d/c. The angular distribu- 
tion is no longer simply given by (]^),(Q), but corrections 
are down by order (d/R) 2 . 

The result (^) is very powerful. The anisotropy a is in- 
dependent of the diffusion coefficient, as long the approxi- 
mation of isotropic source and spatially uniform diffusion 
coefficient is adequate. As a result, the anisotropy is in- 
dependent of the CR energy, and the strength and struc- 
ture of the magnetic field, so long as the propagation is 
diffusive. This has several important implications. First, 
knowledge of a determines the ratio R/Tc where R is 
the distance between source and Earth and T is the time 
between emission and observation at Earth, for instanta- 
neous sources (see below for extended sources). Second, 
UHECRs at all energies can be combined to improve the 
statistics in the determination of a without losing its util- 
ity for determining R/Tc. The highest energy events can 
be excluded, or treated separately, since fields may not be 
strong enough that the highest energy events propagate 
diffusively. An energy dependence of the anisotropy sig- 
nals a breakdown in this simple description, e.g., due to 
large scale structure in the magnetic fields. The angular 
distribution is known to be grossly isotropic, but a quan- 
titative fit to the form (^|) with respect to an arbitrary 
direction is not available. 
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For an AGN, we must integrate eqn (0) over the 
range of possible propagation times. There is a negli- 
gible contribution from times prior to the arrival time 
of the diffusion front, so the minimum propagation time 
is train = Max{i? 2 /(3£>), T ff}, where T a g is the time 
since the AGN turned off its UHECR production. The 
maximum propagation time is t max = MinjrGZif), Ton}, 
where tqzk ~ 50 Myr for 10 20 5 eV and T on is the time 
since the AGN turned on. Idealizing the emission to be 
uniform with a rate uq/t, we have: 

[8ttD(E,B)] 3 / 2 n ax 

since L> - T 1 / 3 for Kolmo gorov diffusion, and shock 
acceleration results in a spectrum at the source of E~ p 
with p slightly greater than 2, the resultant spectral in- 
dex is close to the ~ —2.7 which is observed. The data 
above E ks 10 20 eV is inadequate to decide whether there 
may be a transition to Bohm diffusion with its steeper 
spectrum. The anisotropy is the flux-weighted average of 
R/2tc, so a a R/(6t m inc), introducing some energy de- 
pendence if t m i n is determined by diffusion rather than 
the AGN turnoff. 

These results allow us to constrain the source of UHE- 
CRs observed at Earth. We must demand that the num- 
ber density of UHECRs from the proposed local source be 
equal to or greater than the accumulated number density 
from all UHECR sources in the rest of the Universe: 

n E >T s nH'\ (6) 

where n is the total number of UHECRs produced by 
an average source, T s is the number of sources per unit 
volume and unit time, and w 10 4 Myr is the age of 
the Universe. For a bursting source such as a GRB, tie is 
given by (^), with no = nn to allow for an enhancement 
factor k if the local source is unusually bright; for an 
extended source such as an AGN, he is given by (||). 

We now show that GRBs are unlikely to be respon- 
sible for UHECRs, as follows. The rate of GRBs is 
r s = 0.2 x 10"3 Tsch/b Mpc~ 3 Myr" 1 ||, where 
b = O.Ol&o.oi < 1 is the beaming fraction of GRBs; 60.01 
may be as small as 1/2 @|. To satisfy (§) with (|) 
thus requires DT < 0.0012 (6 . i K/r Sc/l ) 2/3 Mpc 2 , while 
R 2 < 6DT for the diffusion front to have reached Earth, 

implies i?<80 (60.01 ^/^SchY^^pc. Thus to dominate 
the rest of the GRBs in the Universe, a single GRB would 
probably have to be in our galaxy! But a magnetic field 
structure capable of satisfying the isotropy constraint in 
that case could confine cosmic rays below the GZK en- 
ergy from earlier GRBs in the galaxy, and the accumu- 
lation problem is not avoided. The basic problem is that 
the observed rate of GRB's is too high to allow Earth's 
environment to be dominated by the output of a local 
GRB. 

AGN's are much rarer than GRBs and the local dom- 
inance requirement is readily satisfied for the values of 



D and the time scales under consideration here. The 
local flux (||) depends on t~*„ — tmax, which we es- 
timate to be ~ 0.1 - 0.2 Myr~ 1/2 . The AGN rate is 
« Pagntagn/H- 1 , where pagn ~ lO^Mpc^Myr -1 
is the number density of AGNs. (^|) is satisfied with- 
out assuming the UHECR output of the local source is 
stronger than average. It is helpful that taqn cancels 
out, since it is uncertain. 

We now ask for a suitable candidate source. Defining 
D = 0.25 D K Mpc 2 /Myr and taking t min < t rnax <50 

Myr, the UHECR diffusion front is closer than about 
1/2 

<5D^ Mpc. It would be very difficult for the source to 
be significantly farther than this. Thus M87 at 18 Mpc, 
which has been mentioned as a possible single source of 
UHECRs, is too far away if we require diffusive propaga- 
tion. Cen A, at 3.4 Mpc, is by far the nearest powerful 
radio galaxy; see ref. fig] for a comprehensive review. Al- 
ready in 1978, Cavallo |l6| pointed out that the size of 
its radio lobes and the strength of its magnetic fields sat- 
isfy the Hillas criterion (Tt]] for acceleration of UHECRs. 
However at b = 20° and I = 310°, Cen A is in the blind 
direction of modern UHECR detectors which are located 
in the northern hemisphere, and thus can only be the 
source if the UHECRs propagate diffusively. Until re- 
cently there was a prejudice that extragalactic magnetic 
fields are of order nG and deflection of UHECRs is neg- 
ligible, so Cen A was not considered acceptable. 

The observed total luminosity of Cen A is now about 
10 43 erg/s, of which about half is high energy Jl5| |, so 
if e is the efficiency of UHECR production compared 
to photon production extrapolated to 10 19 eV using 
equal power per decade, we estimate E 2 dN/dE/dt as 
3el0 53 eU/s for 10 19 eV UHECRs. Using (§) and i m ^ 2 - 
tmax = O.lMyr ' to be conservative, we obtain the 
energy-weig htcd flux per str at Earth E 3 dN/dE « 
2el0 25 eV 2 m~ 2 s _1 sr _1 , easily consistent with the ob- 
served value of 10 24 5 eV 2 m~ 2 s _1 sr _1 . Furthermore, 
the time scale for evolution of AGNs is w 10 Myr in 
general, and of Cen A in particular fl5| , so that it was 
likely to have been a powerful AGN within the most re- 
cent GZK time. If so, its luminosity would have been 
w 10 44 - 10 45 erg/s and the UHECR production effi- 
ciency need only be of order a percent or less. 

The diffusion coefficient for E = 10 19 eV, B uG = 0.3 
and Am P c = 0.3, is 0.078Mpc 2 /Myr, so the typical 
UHECR arrival time from Cen A is 25 Myr. The path- 
length of a 10 19 eV CR is « 8 Mpc, many times the diffu- 
sion length D/c « 1/4 Mpc, so the diffusive approxima- 
tion applies; the Fly's Eye event barely satisfies the Bohm 
diffusion criterion. The anisotropy a « R/(6t m i n c) is 
predicted to be w 0.07 for 10 19 eV, or less if the AGN 
quit emitting UHECRs more than w 25 Myr ago. 

Having identified Cen A to be a good candidate source 
for UHECRs, we ask which new physics explanations 
are consistent with our three constraints. Lorentz In- 
variancc violation, such that even 3 x 10 20 eV protons 
are below threshold for photo-pion production, is ac- 



4 



ceptable. The Z-burst model involves messenger neu- 
trinos produced throughout the Universe, for instance in 
distant AGNs, but it is functionally a model with dis- 
tributed sources, namely the target dark matter neutri- 
nos, and will be discussed with models of that type below. 
That leaves two GZK-evading messenger models: neutri- 
nos which interact strongly with nuclei in Earth's atmo- 
sphere, or light gluino-containing baryons whose thresh- 
old for photo-pion production is above 3 x 10 20 eV. These 
models are acceptable only if nucleon emission by the 
accelerator is <l// acc times the messenger emission. If 
the messenger neutrinos or hadrons are produced by elec- 
tromagnetically accelerated protons, as is most conven- 
tional astrophysically, this condition demands consider- 
able nucleon attenuation in the source. It is noteworthy 
that this nucleon attenuation requirement implies that, 
if the source were an AGN, it would necessarily display 
the spectral characteristics found to be associated with 
candidate distant sources in the analysis of jl8). If the 
hints of clustering of events and directional identification 
with powerful distant matter-enshrouded radio quasars 
become real signals with better data, these models would 
be favored. If UHECRs are protons from Cen A, these 
hints are mere fluctuations. 

The observed isotropy of "low" energy UHECRs rules 
out all models in which the UHECR sources are propor- 
tional to the galactic DM distribution. The DM distribu- 
tion is fairly well known and there is little flexibility in ad- 
justing its typical clustering scale. AGASA compared the 
anisotropy in the angular distribution of 581 events above 
10 19 eV with that of two popular dark matter distribu- 
tions centered on the Milky Way || and in both cases 
found poor agreement with the predicted anisotropy (re- 
duced x 2 > 10.0). Above 4 x 10 19 eV the statistics are 
inadequate to make this analysis, but with the assump- 
tion of a single UHECR population, that is not necessary 
to exclude such models. Sec Jl9| for a numerical study 
which reaches the same conclusion. 

Relaxing the conditions on the source distribution, as 
may be appropriate for a new population of astrophysical 
sources such as magnetars Q, topological defects, or eV 
neutrinos, does not allow this constraint to be evaded. To 
see this, consider a simple toy model in which all sources 
are on a shell of radius d around the galactic center, with 
d 8.5 kpc, the distance of Earth from the galactic 
center. This source distribution gives a dipole anisotropy 
of 0.34/(d/50kpc). By increasing d the anisotropy can 
be decreased, however large values of d are unnatural; 
e.g., an isothermal halo with a core radius of ~ 1.5 kpc 
has (d) ~ 10 kpc. An upper limit on the anisotropy 
implies a lower limit on d, but increasing d makes the 
local-dominance requirement problematic. 

To conclude, we have argued that Earth at this epoch 
is not exposed to a generic UHECR spectrum. The ab- 
sence of a cutoff at the GZK energy is a reflection of a 
coincidental position in space-time relative to the near- 
est source. The most straightforward explanation which 
survives the analysis here is that most UHECRs reaching 



Earth come from a single, unusually powerful AGN at a 
distance of a few Mpc. Cen A, a powerful radio galaxy 
at 3.4 Mpc is an excellent candidate. This requires dif- 
fusion in a magnetic field, probably in the few-tenth fib 
range, which is consistent with observational and theoret- 
ical constraints. We predict the anisotropy of UHECRs 
with E > 10 19 eV is of order 7% or less. Our analytical 
discussion should be followed up by numerical studies, 
especially to refine our estimates of the magnetic fields 
needed to diffuse the highest energy CRs. 
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